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1. INTRODUCTION 

Here we review recent and new developments which elucidate the symmetry and 
the description of the magnetic and dielectric states of multiferroics using order 
parameters. After some examples where the magnetoelectric (ME) behavior is 
relatively simple, we discuss a proposed generic phase diagram for the so-called 
"125" systems, RMn205, where R is a rare earth. The most important consequence 
of the phenomenological theories we develop is to provide a general framework for 
understanding the magnetic and dielectric properties of these materials and how these 
properties combine to produce the interesting ME phenomena. 

Briefly, this article is organized as follows. In Sec. 2 we discuss the 
characterization of the magnetic structure obtained from symmetry arguments. Here 
we discuss briefly a simplified version of the group theoretical approach (known as 
representation theory) which is supplemented by less well-known arguments involving 
the use of inversion symmetry. As examples we consider Ni3V208 (NVO), TbMnOa, 
and RbFe(Mo04)2 (RFMO) and discuss the introduction of order parameters (OP's) 
to characterize the magnetic symmetry. We then give a brief review of how symmetry 
restricts the form of the ME interaction when it is written in terms of both magnetic 
and dielectric OP's. In Sec. 3 we give a detailed discussion of how these concepts 
enable us to construct a generic phase diagram for the 125 family of multiferroics, 
which does not rely on a knowledge of the details of the microscopic interactions. 
Section 4 contains an application of the theory of critical phenomena to the 125's, and 
the paper is briefly summarized in section 5. 

2. SYMMETRY AND MAGNETIC STRUCTURE 

Here we give a simplified review of the role of symmetry in determining the structure 
of the magnetically ordered phase which develops at a continuous phase transition. 
This subject is of ancient vintage, being discussed about 60 years ago by Landau 
(see [T]). However, some reviews which discuss the analysis of diffraction data[H |3] 
overlook the importance of inversion symmetry in reducing the number of parameters 
needed to describe the ordered magnetic structures. For multiferroics this was first 
corrected quite recently by Lawes et al. [4] , by Kenzelmann et at j5, and in more detail 
by Harris [6], which we follow here. Formal treatments appeared some time ago [7]. 
Recent papers include Schweizer et aL[8] and Radaelli and Chapon[9]. 

We start by assuming that the paramagnetic phase is characterized by a primitive 
unit cell with n^ magnetic sites. The Landau expansion of i^2, the magnetic free energy 
at quadratic order in the spin components, is 

^2 = E E J2^X-\<i)U.;r'J3S^icL,rrS^{ci,T') , (1) 

q t,t'=1 a, (3 

where % is the wave- vector dependent susceptibility matrix and 

5„(R,t)= ^5„(q,T)e^'i-^, (2) 

q 
where Sa (R, t) is the a-component of spin of the rth magnetic site in the unit cell at R 
and Sa (— q, t) — Sa (q, t)* . For each value of the wave vector the inverse susceptibility 
has Srir eigenvalues (which may or may not be distinct from one another). At 
high temperature T all these eigenvalues are positive and the paramagnetic state is 
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thermodynamically stable. As T is reduced through a critical value, T^, one eigenvalue, 
Ac(qc), at some wave vector q^ (and wave vectors equivalent to it by symmetry which 
comprise the star of qc) approaches zero, signaling an instability of the paramagnetic 
phase to the formation of long-range order at the critical wave vector qc associated 
with this critical eigenvalue. The actual value of qc is determined by the microscopic 
interactions. Since these interactions are not well known, we regard the wave vector as 
an experimentally determined parameter. The degeneracy of this critical eigenvalue 
Ac is riqM, where Uq is the number of wave vectors in the star of q and A/" is the 
dimensionality of the irreducible representation (irrep) of the symmetry group (the 
so-called "little group") which leaves the wave vector invariant. (For ferromagnetic 
Ising, x-y, and Heisenberg models Af assumes the values 1, 2, and 3, respectively.) 
To avoid technicalities, in this section we consider the simplest case, TV = 1. This 
case is simple because then we can use the familiar principle that the eigenvectors 
of a matrix (here the inverse susceptibility) are also simultaneously eigenvectors of 
operators (here the symmetry operations Oi of the space group which leave the selected 
wave vector invariant) which commute with each other and with the matrix. In this 
way we avoid using the full apparatus of group theory and the reader need not know 
anything at all about "irreps". We now illustrate this idea and show how inversion 
symmetry introduces further simplifications for three recently studied multiferroic 
magnetic materials, whose lattice structures are shown in figure [T] and whose positions 
(except for RFMO where the Fe ions form a Bravais lattice) are given in table [TJ [j 



Figure 1. (Color online), (a) The six Ni sites in the unit cell of NVO. Circles 
represent "spine" sites and squares represent the "cross-tie" sites. The axis of the 
two-fold rotation about the a-axis is shown. The glide plane is indicated by the 
mirror plane at 2: = | and the arrow above rric indicates that a translation of 6/2 
in the b-direction is involved, (b) The four Mn sites (small circles) and four Tb 
sites (large circles) in the unit cell of TbMnOs. The glide riia is indicated by the 
mirror plane at a; = ^ followed by a translation of fe/2 along the b direction. The 
planes at 2 = 1/4 and z = 3/4 are mirror planes, (c) RFMO, where large balls 
are Fe spins 5/2 on a stacked triangular lattice, small balls are oxygens. Mo ions 
are inside the oxygen tetrahedra, and the Rb ions are not shown. 




(C) 




2.1. NVO 

For NVO the incommensurate (IC) wave vector for magnetic ordering is[TTl [12] 
q = 0.28(27r/a)a. Thus the space group operations Oi which leave the wave vector 
invariant are generated by 2q, a two- fold rotation about the a-axis and passes through 



I We interchangeably denote the a, b, and c axes as x, y, and z, respectively. 
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Table 1. General positions [10^ (given as fractions of lattice constants) within 
the primitive unit cell for NVO (space group Cmca) and TbMnOa (space group 
Pbnm). Here r+ = r + 1/2, 2q is a two-fold rotation (or screw) axis, and nia is 
a mirror (or glide). 



NiaVzOg 


TbMnOa 


Er = {x,y,z) 


2cr = {x,y+,z+) 


Er= {x,y,z) 


2aV= (a;+,y+,z) 


2ar = ix,y,z) 


2f,r = {x,y+,z+) 


2cr = {x,y,z+) 


26r=(x+,y+,z+) 


Jr = {x, y, z) 


m^r = {x,y+,z+) 


Xy ~ (x, y , z) 


niar = {x+,y+,z) 


rUaT = (x, y, z) 


mbV = {x,y+,z+) 


riicV = {x,y,l+) 


nn,r = {x+,y+,z+) 



the origin and rUc, a glide operation which takes c into — c followed by a translation 
through (6/2)b. Thus the critical eigenvector (which is the spatial Fourier transform 
of the spin distribution) must not only be an eigenvector of the inverse susceptibility 
matrix, but it must also simultaneously be an eigenvector of both 2^ and rric. Since 
\^aY — 1? the eigenvalues of 2^ must be A(2a) = A = ±1. Since [rraj^ is a translation 
along the h axis, the eigenvalues of nic must be A(toc) = A' = ±exp(i6(7b/2) = ±1. 
Thus, if we assume continuous transitions, there can only be four distinct symmetries 
of ordered phases, corresponding to independently selecting the eigenvalues of 2^ and 
rUc- The corresponding eigenvectors must be of the form 

S{q, 1) = (ai, a2, aa)^ , S{q, 2) = A(ai, -a2, -as)? , 

S{q, 3) = AA'(-ai, a2, -aa)^^ , S{q, 4) = A'(-ai, -aa, +03)/^^ , 

S{q, 5) = ([1 + A]a4, [1 - A]a5, [1 - AJae) , 

^(g, 6) = - A'([l + A]a4, [1 - A]a5, [1 - A]a6)e' , (3) 

where ^ — eyi^{iqxa/A) and the a„ assume arbitrary complex values. To check this 
note that under 2q sublattices #1 and ^2 are interchanged as are #3 and #4, whereas 
under rUc sublattices #1 and #4 are interchanged as are #2 and #3. Note that 2a 
changes the signs of the b and c-components of spin, while rric changes the signs of the 
a and b components of spin since spin is a pseudo- vector. This type of analysis, known 
as representation theory, is well-known and widely used. However, less well-known and 
often overlooked (as documented in [B]) is the fact that in these multiferroic systems 
the free energy must be invariant under the inversion symmetry T possessed by the 
lattice [7]. One can then showjH [T51 [T^ IB] that this symmetry fixes the phases of the 
a„: for A = A' = 1, apart from an overall complex phase factor, ai and a^ must be 
pure imaginary and a2 and ai must be pure real. For other irreps [i. e. for the three 
other choices of the eigenvalues \{Oi)] one has analogous results. If ([T]) is generalized to 
include terms of fourth order in the spin variables, then a mean-field analysis for T near 
Tc shows that the overall amplitude of the spin wave function varies [proportionally 
to (Tc — T)^"], but the ratios among the a„'s are nearly temperature independent. 
Therefore we replace a„ by cr(q)a„ and require the normalization ^ \o-n\'^ = 1- Thus 
the temperature dependence is incorporated in the order parameter a. If we require 
that a4, say, be real, then the freedom to fix the overall phase is taken into account 
by allowing the order parameter to be complex, as one would expect, since the origin 
of the IC ordering is not fixed, at least within i^2- It should be noted that the order 
parameter inherits the symmetry of the spin functions, so that 

2a<J — \{2a)a — Act , mc<J = X(jnc)(T = A'ct , 2a = a* . (4) 
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In the analysis of diffraction experiments one tries to fit the structure assuming 
in turn each of the four symmetries. In so doing one has not Sn,- = 18 complex- valued 
fitting parameters, but rather the 4 or 5 a's of ^ (depending on which symmetry one 
is considering). However, the use of inversion further reduces the number of fitting 
parameters by half since their phases are fixed [12] . 



Figure 2. Dielectric and magnetic phase diagrams of NVO (left, from |lllHl[T2l ). 
TbMnOs (center, from HUEI), and RFMO (right, from [T71[TS]). In the dielectric 
phase diagram the direction of the spontaneous polarization (if any) is indicated. 
For NVO Tc ~ 4K, T< ^ 6K, and T> ^ 9K. In the magnetic phase diagrams P 
denotes paramagnetic, HTI denotes a dominantly collinear IC phase with a single 
OP, LTI is a dominantly elliptically polarized phase with two OP's, and IC-TRI 
denotes the IC stacking of triangular lattice antiferromagnets. 



Diel. 



Mag. 



Diel. 



NVO T=2 Tc 



LTI HTI P Mag 



lie II ( 



TbMnO T= 



10 30 40 



C LTI HTI 



Diel. 
RFMO 

Mas. 



IC-TRI 



The left phase diagram in figure [5] shows the experimentally observed sequence 
of magnetic phases of NVO. As T is lowered, the first ordered phase to appear 
(at T> « 9K) is the HTI phase, which has a single eigenvector associated with 
predominantly collinear sinusoidally modulated order. Analysis of experimental data 
indicates that the eigenvalues of this phase are[12] X{2x) = —1 and X{mc) = +1. 
At a lower temperature (r< ~ 6K) the LTI phase appears with an additional order 
parameter associated with dominantly transverse spin order and with X{2x) — +1 and 
A(7Tic) = +1. The magnetic free energy which describes the development of these two 
successive ordering transitions is of the form[T21 [H [H] 

Fm = a{T - r>)|aHTi|' + h{T - T<)\cnjTi\'' + 0{n^) , (5) 

where a and b are constants and r> and T< are the respective temperatures (when 
nonquadratic terms are ignored) at which cthti and ctlti become nonzero. The 
unwritten terms in Fm, which are quartic in a, favor fixed spin length. Thus ctlti is out 
of phase relative to cthti and the spins thereby order in a spiral structure[T2 l [T5| . IT6]. 



2.2. ThMnOs. 

The case of TbMnOs is almost identical to that for NVO. Here the IC wave vector is 
of the form fO. q. 0) [roi [5]. so that the symmetry operations which leave it invariant 
are generated by the glide nia and the mirror mc. The eigenvalues of rric {rria) are ±1 
(±A), where A = exp(i7r(7). For A(mc) = 1 and \{ma) = A one has 

S{q, 1) = aii — a2J - ask , S{q, 2) = aii + a2J + ask , 

S{q, 3) = -aii + a2J - a^k , S{q, 4) = -aii - a2J + ask , 

S{q, 5) — a^k , S{q, 6) = —a^k , 

S'((7, 7) = asfc , S'(q, 8) = -a4fc , (6) 

where the a„ are arbitrary complex numbers. As for NVO one can now require that F2 
be invariant under I. In this case the result is that apart from an overall phase factor, 
an for n = 1,2,3 are real, a^ is an arbitrary complex number, and a^ = a^ [5l [131 [6] . 
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When inversion symmetry relates sites within the same WyckofF orbit of the operators 
of the group of the wave vector (as it does for NVO) , the complex phases get fixed, 
whereas when inversion relates sites in different Wyckoff orbits (as for the Tb sites 
here), the previously independent amplitudes of the two orbits are now related. Again, 
we replace a„ by cr(q)Q;„, so that the temperature dependence of the spin function 
is essentially contained in the order parameter cr(q) and the complex phase of the IC 
ordering is incorporated in the arbitrary complex phase of a{q) which transforms as 

macr(q) = A(r7ia)cr(q) , r7iccr(q) = A(TOc)cr(q) , I(j{q) = cr(q)* . (7) 

The center phase diagram in figure [2] shows the sequence of magnetic phases 
of TbMnOs. As the temperature is lowered (through r> = 40K) the first ordered 
phase to appear is the HTI phase in which the single eigenvector associated with 
predominantly coUinear order appears with X(mc) = — exp{iTrq) = —A and X{ma) = 1. 
At a lower temperature (T< Ri 30K) the LTI phase appears with an additional 
order parameter associated with transverse spin order and with X{mc) = —A and 
X{ma) = — 1. The phenomenology of the magnetic phase diagram of TbMnOa is very 
similar to that of NVO. 

2.3. RFMO 

The magnetic Fe ions in RFMO form triangular lattice planes which are stacked 
directly over one another, as shown in figure file) [17] . Below T = 180K but above 
the magnetic ordering temperature (at Tc = 4K) the lattice has P3 symmetry [20], so 
that the only symmetry operation (apart from T) is a three-fold rotation Tt about the 
c-axis, which is perpendicular to the triangular lattice plane. At low fields, the spins 
within a single triangular lattice plane form a 120° structure and as one moves from 
one plane to the next the spins are all rotated through an angle Scj) — QcC, so that the 
component of the IC wave vector along c is Qc [H]. To generate the 120° structure, 
the in-plane component of the wave vector must be chosen to be at the corner, X, 
of the Brillouin zone of the triangular lattice, i.e. q = X + qcC. Then the symmetry 
operations Oi which leave the wave vector invariant are TZ and TZ^^. {TZ takes X into 
a vector equivalent to X.) We thus end up with a one-dimensional irrep F and its 
complex conjugate F*. The spin distribution is given bv[l8l[6j 

S(r) = [o,{q,)(z + zj) + 02{q,)(z - zj)]e''i- + c. C. , (8) 

where /x = exp(27ri/3). The order parameters transform as 

T^t^niQz) = M"'^n('?^) , It^niqz) = (^3-n{qz)* ■ (9) 

The magnetic free energy up to order a'^ is 

F^{T- T,)(T^ + ua^ + v\ai{q,)a2{qz)\'' , (10) 

where ar^ = \cri{qz)\'^ + \cr2{qz)\'^ , and u and v are constants (with u positive). It is 
found [IH] that only one of the two order parameters is nonzero in a single domain, 
from which we deduce that v must be positive. (This conclusion is confirmed by the 
appearance of ferroelectricity, as we will see in a moment.) 

2.4. Magnetoelectric Interaction 

Here we describe the ME interaction which leads to a spontaneous polarization induced 
by magnetic ordering which breaks inversion symmetry. For this purpose we show the 
dielectric phase diagrams of the three systems under consideration in figure [H 
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We write the free energy as 

F = FM+FE + Vint, (11) 

where Fm (Fe) is the magnetic (dielectric) free energy and l^nt is the ME interaction 
which is responsible for the magnetically induced ferroelectricity. 

We first consider NVO[4lfT2] and TbMn03j5il6j. Both have two magnetic ordered 
phases, the high-temperature incommensurate (HTI) phase at higher temperature 
(T> > T > T<), described by a single order parameter cthti for which spins 
are predominantly confined to the easiest direction, and the low-temperature 
incommensurate (LTI) phase (for T < T<) in which a new order parameter ctlti 
appears, describing ordering transverse to that of cthti- The order parameters are out 
of phase (to minimize the fourth order terms in the magnetic free energy) [T^], and 
thus give rise to a magnetic spiral. These order parameters transform as specified by 
(lU and ([7]), respectively. 

We have Fe = {1/2)xe^P'^, where xe is the dielectric susceptibility and P is the 
polarization vector. Because there is no tendency for ferroelectricity to form in the 
absence of magnetic ordering, xe never gets large. In the absence of ME coupling, 
the equilibrium value of P is zero. The ME interaction has to conserve wave vector 
and be invariant under time reversal. At lowest (quadratic) order in cr, it therefore 
must be of the form T^„t ~ <^{<l)<^{~<l)P = crcr*P- In the present situation, the two cr's 
can not both be HTI or LTI, because then Vi„t would not be invariant under spatial 
inversion. So 

V^nt = y^[c-yCrHTl(g)crLTl(g)* + C*crHTl(g)*CTLTl(g)]-P7 , (12) 

7 

and to be invariant under inversion we must have c-y = ir^^ where r-y is real, so that [3] 






^nt = « 7 ^r-y[o-HTl(g)o-LTl(<?)* - crHTl(g)*crLTl(<3')]-P7 (13) 



7 






2sin(0HTi - '/'ltOIo-htiO'ltiI > ,''7-P7 : (14) 



7 



where cthti = IcthtiI exp(i0HTi) and similarly for ctlxi- The transformation 
properties given in ^ and ([7]) for the order parameters under the mirror and glide 
operations then imply that r-y in (|13p is only nonzero for 7 = 6 for NV0[4j [121 HSl E] 
and 7 = c for TbMnO-i [51 [T51 IBj . The fact that P is proportional to IcthtiCTltiI has 
been experimentally verified for NVO[22]- 

For RFMO the argument is slightly different. There ^ indicates that 
o'i{qz)<^2{Qz)* is invariant under inversion (which changes the sign of P). Thus ^ 
implies that the ME interaction quadratic in a, which conserves wave vector, is|18i[6| 






Vint = } r^[\a,{qz)\' ~ \a2iqzr]Py , (15) 






where r-y is real valued. Since the square bracket is invariant under the three-fold 
rotation TZ, P^ must also be invariant under TZ. So at this order r^ can only be 
nonzero for 7 = c, as is observed[TB]. At higher order [5^ a transverse polarization is 
in principle possible. Note that TZ{Px — iPy) = lJ-{Px — iPy) and T^-ctiCTj = /^^ctictj- 
Then one can have an ME interaction of the form 

^it' = c[|ai(q,)p - 1^2 (<Z.)|']^l(<Z.)^2 (<?.)* (P. - iPy) + c.c. . (16) 
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However, the fourth order terms written in (|10p select |(Ti(92)| = |c2(9z)| if ^ is 
negative and ai{qz)(J2{qz) = if u is positive. In either case VJ^^ does not come into 
play. Since the ordered phase is ferroelectric, we deduce that v is positive and that 
only Pc is nonzero. Then, within mean field theory, Pc is proportional to (|cr|"}, as is 
the intensity of the magnetic Bragg peaks. This is experimentally confirmed|18|. |§| 



3. 125's 



Figure 3. (Color online) (a) Two views of the lattice structure of the 125's. 
(b) Symmetry operations of space group Pbam. nia denotes a mirror or glide 
operation, 2^ is a two-fold rotation or screw operation and r+ = r + 1/2. 




(b) 



Er =(x,y,z) 

lY-(x,y,z) 

2^ r ^(x,y,z} 

in^r-(x,y,z) 



2^r =(x+,y+,z) 
2^r ^(x+,y+,i) 
m^ r =(x+,y+,z) 
mj^r =(x+,y+,z) 



We now consider the "125" orthorhombic (space group Pbam) family RMn205 
(RMO), where R=Y, Ho, Er, Dy, Tb, Tm. Their lattice structure and the 
corresponding space group operations are shown in figure [3l The paramagnetic unit 
cell of the RMO's contains 12 potentially magnetic ions: 4 Mn'^+, 4 Mn'*+ and 4 R'^+. 
Experiments show that all the RMO's exhibit magnetic spin density wave ordering, 
with a wave vector q which undergoes a sequence of phase transitions|24j-p9]. To 
discuss these phases we introduce the notation q = {U,0,V)n, which we abbreviate 
as ([/, V)n (in figure [5] these are denoted by UVn)- If C/ = C (C/ = /), then ^ = 1/2 
(g^ = 1/2 - S) and ii V = C {V = I), then q, ^ 1/4 (q^ = 1/4 + e), where the 
wave vector is in reciprocal lattice units and S and e are of order 0.01 and depend on 
temperature. V — X includes the cases when e y^ and when e = 0. The subscript 
n, if it is given, indicates indicates the number (1 or 2) of OP's, see below. As the 
temperature T decreases, all the RMO's (with the possible exception of R=Dy) first 
order below Tc (« 45 K), into an incommensurate (/,/) phase with no ferroelectric 
(FE) order. For YMO (at Tp = 41K)r3Q', 11], ErMO (at Tp = 39K)[32] and TmMO 
(at Tp = 39K)[33]. this paraelectric incommensurate state gives way to an (/, C) phase 
and this phase displays a weak FE moment P along the 6-axis. Below Tc ^ 37 — 39K, 
q locks into a commensurate (CM) value {C,C) and Pb increases significantly [3T]. 
TbMO^ Ml, HoMOlSa |38], and probably DyMO03 go directly from the (/, J) 
phase into the ferroelectric {C,C) phase. At lower temperature (about 10-20K) most 
of the RMO's return to having some kind of incommensurate order [[J. We will not 
be concerned here with these low temperature phases, since their existence probably 
depends sensitively on the details of the spin-spin interactions. As we shall see, the 
behavior of the higher temperature phases can be described by a generic Landau free 
energy. The magneto-dielectric phase diagrams of various 125's are shown in figured] 

§ However, critical fluctuations may imply different exponents for Pc and |crp, see Sec. |4] below. 
II This order may be commensurate but with a large unit cell. 
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Figure 4. ME phase diagrams of ErMO[26], TmMOl28l. YM0[3l l35l . 
HoMO[32l l29l . TbMO[27], and DyMO[39l gO] [31] |32]. We do not indicate 
possible phase changes which have a dielectric signature but only a weak magnetic 
signature and hence may represent a minor spin reorientation. In Sec. 3.4 we 
argue that for 40 < T < 44K DyMO is in an (/, /) phase. 
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Based on the symmetry of the OP's we construct a Laniiau theory for the 
various RMO's, which yields a generic phase diagram, shown in figure [5l which is 
independent of the detailed microscopic interactions[21. Each RMO has particular 
coupling constants which determine the wave vector q. Varying these parameters, Jx 
for Qx and Jz for q^ , changes the value of the optimal q at which magnetic ordering 
occurs. The rest of this section is devoted to an explanation of this phase diagram 
(including the definitions of the various phases) and to a discussion of its consequences. 
This analysis is particularly relevant for the RMO's, because the microscopic theory 
of their multiferroicity is somewhat controversial. Our theory provides a unified 
explanation for the various sequences of phase transitions of the magnetic wave vector, 
and explains why ferroelectricity does or does not occur in the various magnetic phases. 
It also explains the occurrence of two distinct spin structures from neutron diffraction 
studies of the CM phase [35l [25]. This phenomenological theory suggests several new 
experiments and makes a number of predictions, which can be tested experimentally. 



3.1. Magnetic Structure of the {I,I)i and (/, C)i Phases of the 125's 

Given the experimental information, we now analyze the various phases in the order in 
which they arise upon cooling from the P phase. The first phase which is encountered 
is of the (/, /) type. Since qz = 1/4 plays no special symmetry role, it is convenient to 
discuss the (/, /) and the (/, C) phases together. Here, the star of q consists of four 
wave vectors, namely, q± — (±(1/2 — (5), 0, 1/4 + /?) and their negatives. Each wave 
vector is invariant under unity and ruy. This symmetry group has two one-dimensional 
(ID) irreps. Pa and P;,, with complex OP's as amplitudes. By symmetry, all the wave 
vectors of the star must have degenerate eigenvalues of ([T]). Therefore, we introduce 
complex OP's, cr+ = CTa(q^ ) and a~ = a^{ci_) associated with irrep Pa at its wave 
vectors qj?% and similarly for P^. Here q^' and q_' (a^jj and q_ ) are defined to be 
the wave vectors at which the (CTa(q)cra(q)*) ((fTf,(q)crf,(q)*)) susceptibility is maximal 
as T — s- Tea {T -^ Teh). Specific basis functions are given elsewhere |41j. where it is 
also shown that they transform as 

myaM±)^\sa,{c^i^), Xa,(q±) = ^±a.(qL^V, 2ea,(qi^)) = ^^^.(-q^^^), (17) 
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Figure 5. (Color) Left: Schematic 3D phase diagram for q near (1/2,0,1/4). 
The top (red) surface represents the phase boundary between the P and (/,/)i 
phases (where both q^: and q^ are incommensurate). Below the blue surface, 
which is a parabola in Jz (depending only weakly on Jx), one has qx = 1/4, in 
phases (7,C)i and {I,C')2- The green surface represents (/,/)i — > {1,1)2 and 
(/,C)i — > {I,C)2 (the subscripts 1 and 2 denote the number of ID irreps which 
order). Below the orange surface, which is a parabola in Jx (depending only 
weakly on Jz), one has qx = 1/2. Right: A cut at constant qx- The (1,1)2 and 
(^i C)2 phases disappear below the orange surface (as qx — » 1/2), where one has 
a 2D irrep. The dashed and dotted lines are possible trajectories followed by 
specific RMO's as the temperature is varied. 
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where Xa = —At, = eyip^inqx) = rj* and k± — rj^ exp{^2Triqz) . 

As one cools from the P phase, one must enter a phase described by a single irrep. 
Arbitrarily choosing this irrep as Fa, the corresponding free energy is 



p{a) _ 



+ c2k+CT-p + C3[i<y+a-f + c.c.]^4,.,i, (18) 

The coefficients ci, C2 and C3 may differ for F^''\ and we 
When q^ 7^ 1/4, this free energy describes the {I,I)i phase 

+ 1 - \rr~\ 



and analogously for Fh 
assume that Tcb < Tea- 

(the subscript 1 indicates a single irrep). In this phase, we have |ct+| = \a~\ if C2 < 0, 
while only one of a^ or a~ orders if C2 > 0. Replacing (T — Tea) by ^(9o)> where go is 
the wave vector which minimizes r{q), the corresponding minimal free energies in the 
(/, /)i phase are given by Fjj — —r{qo)'^/w, with w = 4ci (4ci + C2) if C2 > (c2 < 0). 
If qz is close to 1/4 then the last {Umklapp) term in (fTS]) can lock q^ to 1/4, via 
a weakly first order transition. Clearly, this term arises only when both a^ and a^ 
order, which would now happen only if C2 — 2|c3| < 0. In this case, one again has 
\'^i\ = Wa\ and Fic — -r(l/4)^/w', with w' = 4ci + C2 - 2|c3|. One would then 
have a first order transition from (/, /)i into (/, C)i when Fjj = Fjc- Since r{q) has 
a minimum at qo, we have r(l/4) « r^qo) + q;(1/4 — qo)'^- Thus, the transition would 
occur when r{qo) + a{l/4: — qo)'^ = r((7o)(u''/^)^ • Remembering that r{qo) = T — Tca, 
we have q^ — 1/4 oc {T^a — T)^/^. Furthermore, since q^ = 1/4 is not a special point, 
we expect qg to be a linear function of J^, hence l/A — q^ tx Jz — Jzc, where Jzc is the 
special value of Jz associated with the transition from the P state into the state with 
qz = 1/4. Thus, the transition from (/, /)i into (/, C)i occurs at T = Tp, with 

Tca-TFO,{Jz-Jzcf, (19) 

as shown in figure [5^. This parabolic relation is a mean- field result. 

Now consider the implications of having qz locked to the value 1/4 in the (/, C) 
phases. From (fT8|) we see that for this locking to occur, both wave vectors q+ and q_ 
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Figure 6. (Color online) Phase diagrams (a) for q^ ~ 1/4, based on I I19I I and 
(b) for Qx ~ 1/2, based on (31) and (32), when the (I,I)i-(I,I)2 phase boundary 
(dashed line) is preempted by locking g^, to q^ = 1/2. Parabolas shown as a 
function of J^ (Jz) are weak functions of J2 (Jx)- (c) r±{qx) for nonzero AJx, 
based on (30). The OP associated with each point is given in the box along with 
the parameters which characterize the wave function, as explained in [41| . In (a) 
and (b) the points M and M' are multicritical points that can only be reached by 
adjusting both the temperature and some additional control parameter. 
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must appear (in the same domain). Then, since we do not allow a direct transition 
from the P phase to the (/, C)i phase (we ignore the unlikely case of a multicritical 
point, at which J^ — Jzc), the question is whether or not a single domain of the 
neighboring (/, /)i phase has two wave vectors. From (IT51) . the condition to have two 
wave vectors is that C2 < 0. An alternate scenario would be that C2 > and the two 
wave vectors do not order simultaneously (in the same domain) in the (/, /)i phase. 
In that case, barring the existence of an as-yet-undetected phase boundary, the two 
wave vectors would have to appear in conjunction with the phase transition between 
the {1,1)1 and (/, C)i phases. For the two wave vectors not to be present in the 
{1,1)1 phase would imply that C2 > 0. Then if C2 — 2|c3| < 0, the two wave vectors 
would appear at the {I,I)i -^ {I,C)i phase transition. It would be interesting to 
experimentally determine (following the logic of Ref. [32]) which scenario actually 
occurs, i. e. whether or not the {I,I)i phase has simultaneous condensation at both 
wave vectors. For this purpose, it would be interesting to perform an experiment 
analogous to that of Ref. [42'. Here, since the ME interaction is present, one could 
use an electric field parallel to one of the wave vectors to manipulate the domains. 



3.2. Magnetoelectric Structure of the {1,1)2 and {I, 0)2 Phases of the 125's 

Similarly to NVO and TbMnOs, the second ID irrep Ft, may order upon further 
cooling. In addition to the 'decoupled' free energies F'^'^^ and _F(''\ the total free energy 
now contains many terms which couple the OP's a^ and a^ . We start by discussing 
the (/, 1)2 phase, where all the wave vector components remain incommensurate. Then 
the quartic terms which couple the two sets of OP's are given by 



ff^=C4[|a>+p4 
where A contains the locking terms. 



C5[\cra(^b 



k-C7+|2 



] + A, (20) 



A = c,[{ataty + {a-a-y + (a+a+*)^ + (a-a-*)1<5 <„, <„ 

+ C7[cri(y-{a+(j'^)* + {<Ja<^a)*(^t'^b\^q^-\,q'^\ 

+ C8[^+T-(a+a-)* + {<^Ub)*^t<'a]6qi^l^qi^l ■ (21) 

Notice that so far we have not assumed that q!^ and q^. are identical. If the 
exchange interactions were isotropic, then the inverse susceptibility would be invariant 



Order Parameters and Phase Diagrams of Multiferroics 12 

under a global rotation of all spin directions. Here, in the generic case, we have small 
anisotropic interactions which break this degeneracy and, in principle, would cause 
these critical wave vectors to be slightly different. In this case, the quartic terms A 
can lock the wave vectors of the two modes into equality, as happens for NV0[T2]. 
The mechanism for this locking is as follows. Assume that, say, aa orders first at 
T = Tea-, and for simplicity we first treat the case with only a single wave vector, so 
that, say, {a'^), but not {(J~), is nonzero. In analogy with what happens for NV0[T2]. 
we assume that ql[° is almost equal to q^ , at which the inverse susceptibility Xb^i'i) 
of (Tf, has its minimum. For Tea > T > Tcb this minimum in Xb (l) i^ positive since 
Ufi has not yet ordered. Now, the quartic terms A give rise to an effective quadratic 
term, V2,cff. Since only (cr^) is nonzero, we have 

^2,cff = c,[{at)^{aty + K+)'K+*)']^q(.,,<,(.) , (22) 

where {X) indicates the thermal average oi X . Even before a^ orders, this term gives 
an additional contribution [beyond (T — Tf-b)] to the inverse susceptibility of a^ , but 
only when q)^' — qy . Since this additional term depends on the relative phase of 
the (Ta's and the ab's, the minimization of this term fixes the phase of cr^, reducing 
its symmetry from that of the XY model (two components of a complex number) to 
that of an Ising model. The minimzation always leads to a negative contribution to 
the inverse susceptibility of a^ . If |(cr^}|^ is sufficiently large, this term can thereby 
shift the minimum in the a^ inverse susceptibility from the wave vector q^ (which it 
would have had when A = 0) into equality with q]^' . Also, the star of the wave vector 
associated with a'^ now contains only the two vectors q^^ and — q+ . This scenario 

applies if the wave vectors for aa and ab are close enough to be locked to q^-^ by the 
term V2,eff before reaching the temperature Tcb at which ab condenses. This, in turn, 
relies on the smallness of the anisotropic terms which cause q^ to differ from q^ • 

If both cra(q+ ) and cra(q_ ) condense at T = Tea, then we need to consider all 
the terms in pT|) . In the (/, /)i phase, both ((T+) — xe"^ and (<t~) = xe^^ break the 
symmetry and have well defined phases </> and x (x is a real number). Substituting 
these values into (j2ip the yields a quadratic form in the four real and imaginary 
parts of e^'f'a'^ and e^^a^ , with eigenvalues 2x^[cq ± (07 + Cg)] and 2x^[cq ± (cy — cg)]. 
Since only one of these eigenvalues is lowest, only one combination of the four OP 
components of a^ orders, and thus we still have an Ising-like ordering into (/, 0)2- In 
any case, we henceforth assume that both OP's have the same critical wave vectors. 

Experimentally, it seems that the phase {1,1)2 has never been observed. Instead, 
the phase with two OP's below (/, /)i seems to be of the {I,C) kind. Therefore, 
we now consider the possible locking of Qz to 1/4, which would correspond to the 
appearance of the (/, 0)2 phase. When g^ is close to 1/4, ([20|) must include additional 
Umklapp terms, which are also consistent with the symmetry of ([TT]) and which lock 
Qz to 1/4. For Qx 7^ 1/2, these are 

Uab = {cga+a+a-a- + cw[{a+a^y + a-a+f] + c.c.}Siq,,i, (23) 

where cg and cio are real. The locking is stronger when two irreps, rather than a 
single irrep as in (fT8|) . are present, because then the additional terms of (j23|l come 
into play. However, in either case, note that this locking requires the presence of both 
wave vectors q+ and q_ . 

Finally, we discuss the ME interactions in the (/, /) and (/, C) phases. In analogy 
with (|13p . the lowest order ME interaction which is invariant under the operations of 
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(HZD is 

^int = irPy ^[fTa(q±)crb(q±)* - CTa(q±) Vb(q±)] . (24) 

± 

Thus, in the (7, J) and in the (I, C) phases, at this order, ferroelectricity requires 
the presence of two order parameters which are not in phase with one another. 
At fourth order in the magnetic order parameters, the ME interaction can lead to 
small spontaneous polarizations in the other coordinate directions, but due to space 
limitations we refer the reader to [3T]. 

3.3. Magnetoelectric Structure of the (C,X) Phases 

This case includes both X = I {qz = 1/4 + e) and X = C {qz = 1/4). Because q is on 
the Brillouin zone boundary [qx = 1/2), the wave vector is invariant under rua and mb, 
and the star of q consists of q and — q. These operations lead to a two-dimensional 
irrep[131[5] and we choose the basis functions as in table XVI of [5]. The actual wave 
function is a linear combination of the two basis functions with complex amplitudes 
o-i(q) and (T2(q). These are the OP's which characterize the magnetic structure and 
they transform as[B] 

mx(Tn{ci) = CnO-„(q) , mycr„(q) = CnCT3-„(q) , 2:cr„(q) = cr3-n(q)* , (25) 

where Cn = (—1)"''"^. Consistent with these symmetries the magnetic free energy up 
to quartic order in a is 

Fm = (t - Tc)[WM)? + k2(q)P] + u[\aM)? + k2(q)P]' 

+ w|CTi(q)cr2(q)|^ + w[cri(q)cr2(q)* + (^2{^)'^l{^)*? 

+ [x{cJi{cif + CT2(q)^) + 2/CTi(q)V2(q)2 + c. c] V-4 > (26) 

where x and y are real. Under the terms quadratic in a and those scaled by u, all 
directions in the four dimensional space of (Ti = (Ti(q) and 02 = 02(0^ are equally 
unstable relative to ordering. However, for qz 7^ 1/4, the fourth order terms select 
|<7i| = |o'2| for ui + 4u < if w is negative, ui — ±jcr2 for tu < if w is positive, 
and criCT2 — otherwise. For qz — 1/4 the terms in x and y are difficult to analyze 
analytically, but in many cases we find that the phases of ui and of CT2 can be chosen 
so that Fm still has minima when either |cri| — |ct2| or cricr2 = 0. 

Now we consider the dielectric properties. At quadratic order in a, since 
Xctictj = ficrj, (fT5|) also applies to the 125's when q = (1/2,0, g^) and then ([25|) 
indicates that r^ is only nonzero for 7 = 6. Including terms of higher order in a [23] 
the ME interaction for the 125's is of the form 

^„t=r,[|aip-|f72p]n + *5K[l'^il'-|'^2p][aia2*-a*a2]F^ , (27) 

7 
where, according to (|25p . the real coefficient r' is only nonzero for 7 = 0. However, as 
mentioned above, (|26p probably allows only either \ai\ = |cr2| or ai(72 — 0, in which 
case the last term in ([?f)) is inoperative. On the other hand, if cri(T2 — (so that, 
say, (72 — 0) and if one applies an electric field, Ea, in the a direction, which induces 
a nonzero value of Pa, the second term in (|27p will induce a nonzero out-of-phase 
value in the order parameter, CT2, that was zero for Ea — 0. Then with (cti) ^ 0, this 
effective linear coupling between Pa and CT2 gives rise to electromagnons[44 j l45 l l46 ] . 

The ME coupling can induce lattice displacements at wave vectors which are 
even integer multiples of the magnetic wave vector [47l |48]. Since the results are 
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particularly simple for the {C,C) phase, where q = (1/2,0, 1/4), we now discuss the 
lowest order interaction in that case. So far we considered a trilinear spin-phonon 
coupling involving cr(q)cr(q)*, which conserves wave vector and therefore couples to 
a uniform polarization. We now generalize this analysis, and consider terms of the 
form (T(q)^ or cr(q)*^, which couple to phonon modes with wave vector ±2q. Within a 
reciprocal lattice vector, this phonon wave vector is equal to the antiferroelectric wave 
vector (0,0,1/2). To construct this interaction we need the site symmetry analysis 
for this wave vector, which is the same as for the wave vector (0, 0, 0) as given in 
Table I of [IS]. There it is indicated that there are 15 63^ (x-like) modes, 15 B2„ 
(y-like) modes, and 9 Bi„ (z-like) phonon modes. An x-like mode, for instance, need 
not involve displacements along the x-axis' rather such a mode need only transform 
like X under the space group operations. Ijl Accordingly, let uaf(7,''") denote such 
a phonon, where 7 labels the symmetry {x, y or z, since we are only interested in 
vector-like modes which carry a polarization) and the index t labels the occurrence. 
We use the transformation properties of ([^5)1 with rriz = Jruxfriy, so that mzUn = cr^- 
Thus the combination (u^ + (t|) is even under rrix and ruy, so that the spin-phonon 
interaction contains the term 

Vsp-ph,z ^ ^[irricrj + crl) + c.c.]uAFiz,T) , (28) 

r 

where r-r is real, so that the square bracket is odd under rn^. Similarly cri(T2 is odd 
under m^ and rUy, so it cannot couple to a vector. Finally {af — cr^) is even under rrix 
and odd under my and it gives rise to an ME interaction of the form 



^sp-ph,y = y'[r:^(CTi - 0-2) +C.C.]uAF{y,T) , (29) 



where r^ is real and we noted that the square bracket is even under rn^ . In summary, 
at this order one can have antiferroelectricity with polarization along either y or z. 

We next analyze the tongue associated with q^ = 1/2. Note that for a critical 
value, Jxc, of the control parameter J^, the two branches [denoted r±(qx, Jx)] of the 
quadratic coefficients r{qx) of the inverse susceptibility are degenerate and are minimal 
at Qx = 1/2, so that r±{qx, Jxc) = '^(0) + a{qx — 1/2)^ + 0{qx — 1/2)**, where a is a 
positive constant. As Jx is varied away from Jxc, a term in r±{qx) which is linear in 
kx = (1/2 — qx) is allowed and generically is of order AJx = Jx — Jxc^^lt- □ The 
symmetry operation mx dictates that the spectrum of the two branches r±{qx) should 
be independent of the sign of kx , as shown in figure [6t, so that 

r±ikx, Jx) = r(0) + akl ± 6fc,(J, - Jxc) ■ (30) 

and for concreteness we assume that the constant b is negative and that Jx > Jxc- 
Symmetry thus implies the existence of equivalent minima at kx = ^b{Jx — Jxc) / i'^o-) = 
k^. Thus at its minimum r±{kx) assumes the value r{kx = 0) — a' {Jx ~ Jxc)"^ , where a' 
is a constant. Accordingly, we can adopt the argument leading to (fT9|) . to the present 
case and obtain 

Tea -Tc CC (Jx - Jxc)"^ , (31) 

^ As explained in |50l . the largest polarization will come from Ta-like modes, which have 
displacements in the ra-direction. 

+ To see the existence of such a term consider the approximation in which, for a system with isotropic 
exchange interactions J(r, r') between spins at r and r', one has, for t ^ t' and q = (1/2 — 5, 0, qz) 
that x~ . / (q) = J(T",a-|-x') exp[7ri — 27ri<5)] + J(x, — a-|-x') exp[— 7ri + 27ri5)]. As long as the sites 
do not sit at a center of inversion symmetry, these two terms will have different amplitudes and will 
give an imaginary contribution which is linear in both 5 and J. 
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where Tc is the phase boundary between the {I,X) and (C,X) phases. (This phase 
boundary is the sohd Hne in Fig. [5)d.) 

The structure of ([50)1 also allows us to discuss the phase boundary Ti^2 between 
the {I,I)i and {1,1)2 phases. For that purpose we compare (|5D|) with P^ (and with 
its analog for i^^''^) and identify r+ with T — Tea and r_ with T — Tct- We thereby 
find that 

-l^c.a ^c.b^ ^0"- yJx ~ •Jxc) '~^ Cy^^x ~ •Jxc) j I'^^j 

where c is a constant. Thus Ti^2 is proportional to (Jx ~ Jxc)^- Depending on the 
parameters, this parabolic tongue can be either narrower or wider than that considered 
above for locking q^ to q^ = 1/2. In the figure we show the former case, since the 
(/, 1)2 phase has not been observed for any of the 125's. 

3.4- Generic Phase Diagram for RMn205 

We now explain how the generic phase diagram of figure [5] describes the various 
RMO's. Since qz = 1/4 is not a high symmetry point, we can not condense from 
the P phase into qz = 1/4 unless we adjust the J's appropriately to reach this higher 
order multicritical point. Since we reject this accidental possibility, the first ordered 
phase we encounter has qz 7^ 1/4. Although q^ = 1/2 is a special value (characteristic 
of antiferromagnetically doubling the size of the unit cell) , the result shown in figure 
\E]p indicates that a continuous transition from the P phase into a (C, /) phase is not 
allowed because it would also involve a multicritical point. For the RMO's (except 
R=Dy which we discuss separately) , experiment shows that the first ordered phase is 
(/, /) and this case is shown in figure [Sj From now on we arbitrarily set Tc^a > Tc^ 
(since we reject the possibility of accidental equality). Consequently we identify that 
the transition from the P phase is into an ordered phase {I,I)i with a single OP 
a^ (except for the star of q). For a single OP, ([M]) provides a phenomenological 
explanation for why this phase is not ferroelectric. As discussed above, we assume that 
in the {I,I)i the phases q^ and qj_ become locked into equality without crossing a 
phase boundary. For the phases with q^ ^ 1/4, experiments have not yet indicated 
whether the two wave vectors qj- occur in separate domains, or whether the true 
state is the superposition, within a single domain, of the two wave vectors. As T is 
further reduced through the (/, /)i phase, a second continuous transition could occur, 
producing a phase (/, 1)2 in which both OP's aa and at, are nonzero (as in NV0[11[5] 
or TbMn03[5j). 

The above description applies for Jz relatively far away from Jzc, i.e. qz relatively 
far away from 1/4. If qz — 1/4, one goes directly from the P phase into the (/, C)i 
phase, which is similar to the (/, /)i phase. Upon cooling, the OP related to the other 
ID irrep tends to order, and one has a transition into the (/, 0)2 phase. This transition 
happens at a higher temperature than that for (/, /)i — > {1,1)2, due to Umklapp 
terms like (|23|) . which enhance the tendency of crb{qx, 0, ±1/4) to order (compared to 
(j{qx,^, qz) with an IC qz). If qz is close to 1/4, one first goes from the P phase into the 
(/, /)i phase, but then the Umklapp terms cause a transition into the (/, C)i phase, 
and one ends up with the phase diagram shown on the RHS of figure [5] 

As the temperature is lowered, each individual RMO follows some trajectory in 
the parameter space. The RHS plot in figure [5] shows possible projections of such 
trajectories. The trajectories, as well as the optimal wave vectors, are assumed 
to have some temperature dependence, which can originate from the elimination 
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Table 2. Values of qx, qz, and n = 1 1/2 — (Ja; | — 1 1/4 — Qz | for T near Tc for various 
RMO's. Positive IT favors locking qx to the value 1/4 in preference to locking q^ 
to the value 1/2. 
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of secondary degrees of freedom, which generate effective temperature-dependent 
exchange coefficients. Note that the whole diagram corresponds to the close vicinity 
of q = (1/2, 0, 1/4), so that this temperature dependence is relatively weak. As shown 
in figure m the real RMO's go directly from the (/,/)i phase into either an (/, C) 
(for R=Er, Tm, Y) or into the (C, C) phase (for R=Ho, Dy, Tb). In the former case, 
we now argue that this phase must be the (/, 0)2 phase: Since the experimentally 
observed phase is ferroelectric, it follows that there must exist two OP's, aa and 
(Tf). Once both order parameters exist, this phase could be either {1,1)2 or (/, (7)2- 
Since the experiments find that qz = 1/4, this must be {I, 0)2- Indeed, we conclude 
that the trajectories for R=Er, Tm, Y are represented by the dashed lines with long 
dashes in the phase diagram. As the same lines indicate, one would then go into 
the (C, C) phase, as indeed observed. At lower temperatures, the trajectories could 
leave the {C,C) phase to the other side of the parabolic 'tongue', and enter a less 
commensurate phase, which could be paraelectric [(/, /)i] or ferroelectric [(1,1)2 or 

(c,/)]. _ 

As indicated by the dashed line with short dashes in the same figure, one can 
also go directly from (/, /)i into {C,C). This trajectory thus describes the RMO's 
with R=Ho, Dy, Tb. In the (C, C) phase, which is ferroelectric, ([77]) indicates that 
|ci| 7^ |cr2|. The quartic term of (pS)) implies that either \ai\ = |ct2|, or one of them is 
zero, so that aia2 = 0. Thus only the first term in (|27p survives and it explains the 
observation [36l [37] that the spontaneous polarization lies along the b axis. Finally, we 
should mention that the fact that different R's follow slightly different trajectories is 
reasonable from the following qualitative point of view. For Tm, Er, and Y the value 
of qx (listed in table [2|) is much closer to 1/4 and therefore is more likely to be locked 
to qx — 1/4 than is that of Ho and Tb. 



For DyMO, experiments have not definitively determined the sequence of phase 
transitions in the wave vector, because the large incoherent neutron cross section 
of the Dy nucleus causes experimental problems. A recent X-ray experiment [40j has 
confirmed the existence[39] of the (C, C) state. The specific heat [37] provides evidence 
that there is a single intermediate phase between this state and the paramagnetic state. 
As argued in connection with figure 5, this intermediate phase has to be an (/, /)i 
phase, because we do not allow the possibility of accidentally hitting the multicritical 
point where the P phase meets the (/, C) (in figure 5a) or (C, /) phase (in figure 5b). 
This proposed phase exhibits a single OP, which is also consistent with the fact that 
DyMO is paraelectric for T > 40K (see figure 4). 

We now return to the phase diagram of figure [S] All the RMO's have q close to 
(1/2,0,1/4) (see table [2]), so they leave the P phase near the apex of the tongue of 
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figure [S^ or[SjD. The effects of a magnetic fieid are explained as follows: it generates 
magnetic moments on the R ions (even above their ordering temperature). Since these 
ions couple to the Mn ions, their moment changes the effective Mn-Mn interactions, 
thus changing the 'control parameters' and the optimal q. This often moves the 
material towards the (C, C) tongue, resulting in a transition from (/, C) ((/, /) when 
paraelectric) back into the CM phase [32l [29]. Pressure[52] has similar effects. 



3.5. Spin structures in the {C,C) phase 

The introduction of OP's leads to a natural interpretation of neutron scattering results 
for the {C^C) phase in YMO. Figure [7] shows the Mn'^+ a-b plane spin components 
in the CM phase of YMO, from the neutron diffraction results of [35, |j and |25j . 
These two structures are obviously similar, and one might ask what symmetry (if 
any) relates them. (This degeneracy was also found in the first-principles calculation 
of [53].) We now show that these two structures are indeed equivalent. 21 To identify 
the symmetry element that relates them note that the structure on the left is even 
under the glide operation nix^ while that on the right is odd under rux. (Here one 
should note that spin, being a pseudovector, transforms with an additional minus sign 
under a mirror operation.) Then (|25p indicates that the structure on the left has 
a2 = 0, whereas that on the right has Gi = 0. Going between these two structures 
corresponds to a rotation in OP space. This equivalence is easily understood when 
OP's are introduced, as done here. Since either ui — Q ov a2 = 0, we conclude from 
the discussion below (|26p . that w + 2v — 2\v\ is positive and both OP's can not order 
simultaneously [6] [21]. This conclusion supports that reached above, namely that since 
the CM phase is ferroelectric, the fourth order terms in (|26|) must select crifT2 = 0. 



Figure 7. (Color online) Schematic diagram of the a and b components of the 
Mn'^^ spins in a single a-b plane of YMO for the CM phase. The glide ma: consists 
of a mirror plane M at x = a/4 followed by a translation b/2 along y. Left: the 
structure given in table III of |25) (with the c-components not shown). Right: the 
structure given in figure 2 of 1351 (who reported zero c-componcnts of spin.) 



M 



YMO M 



j^ ^ ^ ^ -^--^ 



-^ ^ ^ '"^ ^ ^ 



-*-x 3 



To make this identification more quantitative, we consider the magnetic structure 
which H. Kimura et a/. [25] deduced from their neutron diffraction study, which we 
summarize in table [3] Their structure determination was based on an unrestricted 
fit, in which no particular symmetry was assumed. In contrast, our analysis based on 
representation theory assumes that the magnetic structure is characterized by the two 
complex- valued order parameters a\ and (72, with corresponding spin wave functions 



The top (bottom) panel of figure 2 in this paper should be labeled 24. 7K (1.9K). 
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Table 3. Magnetic structures of YMn205 at T = 25K showing the spins vectors 
of the eight Mn''+ and the eight Mn*+ sites within the cell ao X feo X 2co. The 
complete magnetic unit cell is found by antiferromagnetically doubling the cell in 
both the a and c directions. 
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M, 


M^ 


My 


M, 
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M, 


M^ 
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H. Kimura et 


al. 




This work 






H. 


Kimura et 


al. 


1 


This work 




-2.02 


-0.41 


-0.71 


-2.04 


-0.38 


-0.67 


1 


1.74 


0.51 


0.28 


1.72 


0.55 


0.30 


2.20 


-0.40 


-0.24 


2.18 


-0.41 


-0.26 


2 


1.69 


-0.59 


-0.32 


1.72 


-0.55 


-0.30 


-2.06 


0.35 


0.63 


-2.04 


0.38 


0.67 


3 


0.98 


-0.33 


-0.69 


0.99 


-0.30 


-0.63 


-2.15 


-0.42 


-0.29 


-2.18 


-0.41 


-0.26 


4 


1.00 


0.27 


0.57 


0.99 


0.30 


0.63 


2.82 


0.58 


-0.51 


2.85 


0.53 


-0.48 


5 


-1.65 


-0.46 


0.51 


-1.63 


-0.50 


0.56 


-3.07 


0.55 


-0.18 


-3.04 


0.57 


-0.19 


6 


-1.61 


0.55 


-0.62 


-1.63 


0.50 


-0.56 


2.87 


-0.49 


0.45 


2.85 


-0.53 


0.48 


7 


-2.12 


0.74 


-0.09 


-2.15 


0.68 


-0.10 


2.99 


0.59 


-0.21 


3.04 


0.57 


-0.19 


8 


-2.18 


-0.63 


0.09 


-2.15 


-0.68 


0.10 



which are given in [6], but more conveniently in table IX of [4T]. Since we expect 
that cricr2 = 0, our theory would imply that the spin structure should be fitted with 
only one OP component. Indeed, we find that Kimura et al.^s data can be fitted with 
(72 = 0. Optimizing the parameters of table IX of [H] so as to reproduce the spin 
structure of Kimura et al, we found the optimal structure constants to be 

ri = (-0.387, -0.072, 0.09H) , rz = (0.413, 0.078, 0.036i) 

z = (0.257 + 0.049i, -0.081 - 0.017i, 0.031 - 0.063i) . (33) 

With the normalization 2|rip + 2|r2p + 4|zp = 1, the complex-order parameter was 
found to be 



CTi = 5.2698 + i7.3691 



(34) 



(This complex phase can not be explained by a low order anisotropy in the complex cri 
plane.) From table[3]one sees that the structure assuming the validity of representation 
theory is quite close to that of the unrestricted fit of Kimura et al.. The difference 
between these two structures is that our version respects the symmetry one would 
attribute to a structure having only ai nonzero. Thus, in our structure the magnetic 
sublattices are related in pairs, whereas in the structure of [5S] these sublattices are 
almost, but not exactly, related. To characterize the difference between these two 
structures, note that |cri| w 9.1 gives the square root of the sum of the squares 
of the spin amplitudes within the cell of table [H The analogous quantity for the 
difference vector between the two structures is 0.23, indicating that the difference, if 
real, corresponds to an additional order parameter whose magnitude is about 2.5% 
of cTi. As we explained, near the high-temperature limit of this phase one can only 
have either |cri(T2| = or \ai\ — \a2\. Thus, if |cr2| 7^ then we would expect it 
to be of the same order as |cri|. Thus, it seems unlikely that if such an additional 
order parameter would emerge, it would be so small deep in the CM phase, where 
the data were taken. Accordingly, we propose that the actual magnetic structure in 
the (C, C) phase of YMO corresponds to a single order parameter cti. We have also 
identified that the data from [25] on IIoMn205 exhibit the same symmetry: namely the 
(C, C) phase is characterized by the single order parameter cti . Similarly, we identify 
that the magnetic structure of the Mn spins in ErMn205, as reported in [25], is also 
consistent with the symmetry associated with the single order parameter ax. However, 
the phases (px of the x-components of the Er magnetic moments (0.87r and — 0.37r) do 
not agree with the values (tt or 0) corresponding to txi. It would be interesting to 
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check the sensitivity of the data to variation of these phases. It is interesting that the 
structures of all the 125's determined in [53] have 0-2 = 0, even though the structure 
with CTi = represents an equivalent way that magnetic ordering can break symmetry. 
Apparently, the sample preparation (which might create some uniaxial strain) or some 
other experimental detail (stray electric fields?) chooses the structure with a2 — 
in these experiments. It would be interesting to study the cause for this apparent 
symmetry breaking. 

The selection of which OP is nonzero in the (C, C) phase is a result of broken 
symmetry. An electric field along b would order Pf,, and then P^ would select either 
a^ or CTjj depending on the sign of the field. Therefore we suggest cooling the sample 
into the FE phase in the presence of a small electric field along b. Depending on the 
sign of the electric field one should get either the left- or the right-hand panel of figure 
[T] This was indeed confirmed experimentally [55j . (A similar experiment was recently 
performed in TbMnOs^SB ). 

4. Critical phenomena 

All the quantitative results presented above were based on the Landau expansion and 
on mean field theory. Although these theories usually give reasonable predictions far 
away from critical points, fluctuations must be included in the critical regimes. We 
start with NVO and TbMnOa. In these materials, one first goes from the P phase 
into the HTI phase, which is represented by a single complex OP cthti- Since the free 
energy only involves JcthtiP, it does not depend on the phase of this complex number, 
and therefore this transition belongs to the universality class of the XY model, with 
the critical exponents of an isotropic (n =2)-component spin model. The transition 
from the HTI phase into the LTI phase, at r<, is also continuous. A priori, ctlti is 
also a complex number, which would be described by an XY model. However, as we 
discussed after ((2T|) . terms hke [(cthti^lti)^ -|-c.c.] would lock the wave vectors of the 
two order parameters to each other, even before one reaches T< . This lock- in is indeed 
observed experimentally in the LTI phases of NVO [12] and TbMnOs [5]. 

Technically, near r< we have a finite order parameter (cTjjrpj) = ae*". Writing 
also Ctlti = e^'"(5 + ic), the above locking term thus becomes a^(6^ — c^)- Therefore, 
the real order parameters b and c now have different quadratic terms, and only one of 
them (depending on the sign of the overall coefhcient) orders at a temperature slightly 
above the 'bare' T<. As stated above, the fixed length constraint prefers CTjj^j and 
CTlxi to have different phases, which implies that c orders first, and the phases of the 
two order parameters differ by 7r/2. This then yields a helical structure in the LTI 
phase [15l [TBI [12] . Furthermore, this phase relation is also confirmed by the existence 
of a ferroelectric moment in the LTI phase, which would not exist if (pmi = (J)-lti 
(namely if b were to order, rather than c), see ([T4|) . Thus, the transition from HTI 
to LTI belongs to the Ising (n = 1) universality class. Further away from the critical 
point the critical exponents may approach their mean field values 7=1 and f3 = 1/2. 

We next consider the ME interaction, ([T2)l and ([T4|) . Assuming that indeed only 
c orders, we find that near the HTI^LTI transition one can replace ([T4|) by 



Vint = 2rbacPb. (35) 

This immediately implies that the actual order parameter at this transition is not just 
c, but rather a linear combination of c and Pb [54|, [6] . This implies that the dielectric 
constant should diverge near r<, as eb ^ \T — T<^\~'^ , with the Ising susceptibility 
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exponent 7. However, as noted before (J12p . Xe is niuch larger than \T — T^\, and 
therefore the amphtude of this divergent term (related to the amplitude of Pb in 
the mixed OP) can be quite small. It would be useful to search for this divergence 
experimentally. Similarly, we expect that both c and Pb grow below T< as (T< — T)^ ^ 
with the Ising order parameter exponent [3. 

We next turn to RFMO. As discussed in Sec. 12.31 the ordered phase has two 
complex components of the magnetic OP, cri and a2 , and therefore altogether we have 
n = 4 OP components, as described by (jlOp . In fact, this free energy can be written 
as 

F = (T - T,)(|aife)p + |a2fe)n + «(kife)|4 + |a2fe)h 

+ iWiiq.ma^iq.)? . (36) 

This can be viewed as the free energy of two XY models (with OP's Ci and ^2), 
which are coupled by the last term. In terms of the renormalization group (RG), 
this model has two competing fixed points: the isotropic (n = 4) one with v = 2m, 
and the decoupled one with v = OjSj. It turns out that v is slightly relevant near 
the isotropic fixed point, and v is slightly irrelevant near the decoupled fixed point, 
so that as T approaches T< one could follow two scenarios. If u = w — 2m < 0, 
iteration would make it more negative, and one could end up with a crossover from 
the isotropic {n — 4) critical behavior to the asymptotic behavior of two decoupled 
XY models. However, this crossover is very slow. Therefore, one might either observe 
effective exponents close to those of the isotropic (n = 4) critical behavior, or one 
might encounter relatively large corrections to the decoupled critical behavior, due to 
the irrelevant parameter w, which would be renormalized into u(T< — T)~", where a 
is the specific heat exponent of the XY model. Alternatively, if w > then v would 
grow larger under iterations, and one would never reach the vicinity of the stable fixed 
point at -D = 2u + u = 0. In this case, one probably ends up with a slow crossover to 
a weak first order transition. 

The ME interaction in RFMO is given in US]). Thus, Pc - (|cri(<7^)P - \<T2{qz)\^). 
The RHS of this relation represents an order parameter anisotropy. Near the isotropic 
fixed point, this average scales as 

Pc- {\a,{q.)? ^\o2{q.)?) - {\a,\')\ (37) 

where the exponent A > 1 is associated with the scaling of quadratic anisotropy terms 
near the isotropic n = 4 fixed point [58j . However, for this result to hold we must 
have i7ii72 = 0, which arises only if u > 0. As explained above, in this case we expect 
a crossover to a weak first order transition. Thus, as T is increased towards r< we 
would expect a gradual variation from the mean field result, Pc ~ (jcrip), via the 
critical behavior of (|37p . to a weak first order transition. The mean field behavior, 
with A = 1, implies that the FE moment is proportional to the intensity of Bragg 
peaks, as apparently found experimentally [TB^. It would be interesting to check this 
relation close to r< . 

Finally we turn to RMO. As stated, the ordering below the P phase is into the 
{1,1)1 phase, which corresponds to a single irrep, say r^. As seen from ^TS\\ . this 
ordering involves the two complex OP's (T+ and a~ , and therefore belongs to some 
n = 4 universality class. In the (/, /)i phase, where qz 7^ 1/4, the quartic terms in the 
free energy include only those with the coefficients ci and C2. Clearly, this free energy 
is equivalent to the one discussed above for RFMO, yielding only one wave vector if 
C2 > and two wave vectors if C2 < 0. In the former case one probably flows under 
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the RG towards a weak first order transition, while in the latter case one would flow 
towards the stable decoupled fixed point. Thus, the question whether one or two wave 
vectors order is directly related to the nature of the critical behavior. 

The situation changes in the (/,C)i phase, where one also needs to include the 
Umklapp term with C3. Near the decoupled fixed point, this term involves products 
of anisotropies in each of the XY models, and thus it can be shown to be relevant 
[57j . As far as we know, this free energy has no stable fixed point, and one would 
eventually end up with a weak first order transition. However, in the vicinity of the 
isotropic fixed point, where C2 and C3 are small, one could still observe the critical 
exponents of the isotropic n = 4 universality class. In any case, in the generic case 
the phase (/, C)i is reached from the phase {I,I)i via a first order transition, so that 
the critical behavior of the former can only be expected near the multicritical point 
where Jz = Jzc- 

We next discuss the transition into the (so far unobserved) {1,1)2 phase. We 
start with the simple case, where only a^ orders in the (/, /)i phase. As explained 
after ([22)1 , the locking of the wave vectors of (T+ and a'^ fixes the phase of a^ , so that 
the transition into the {1,1)2 phase now involves an Ising-like order parameter. The 
situation now becomes exactly the same as in (|35|) : the dielectric constant e?, would 
diverge with the Ising exponent 7, and Pb would grow in the (/, 1)2 phase with the 
Ising exponent /?. 

The transition from (/, /)i into (/, 0)2 is also weakly first order, since it involves 
a lock- in of q^. However, if the discontinuity is small (as seems to be the case 
experimentally), we can still discuss criticality of the OP's associated with Tt- As 
discussed after (|22|) . this ordering should also belong to the Ising universality class: 
before one reaches this transition one should see et ^ \T — T^i,\~^ and Pb ~ (T^f, — T)^ , 
with Ising exponents. Since A now introduces several additional quadratic terms in 
the fTb's, this transition is expected to occur at a temperature T^^ higher than Tct, 
where one would have the {I,I)i -^ {1,1)2 continuous transition. 

Near the P— > {I,I)i transition (which occurs at Tci), a leading fluctuation 
expansion yields Ae oc {P'^) ex |(CT^)(a'b)|- Since only a^ becomes critical there, 
we expect singularities in e which behave as the energy {\T — Tcil^"") and as the 
square of the OP {{Tci — T)'^^), but with the appropriate effective n = 4 exponents. 
Indeed, experiments [52] show a break in slope at Tci, apparently confirming this 
prediction. This behavior is also expected for other multiferroics and indeed this 
may explain the anomaly seen in the dielectric constant of NVO shown in Fig. 4b 
of [22 . In addition, this anomaly in the zero frequency dielectric constant reflects 
the emergence of a resonance in the frequency-dependent dielectric constant due to 
electromagnons[44l |45l [46] . 

5. Summary 

We have developed a phase diagram to explain the multiferroic behavior of the 
family of 125's systems and have proposed several experiments to explore the 
unusual symmetries of these systems. In view of our current understanding it 
seems unnecessary to invoke the alternate route to multiferroicity proposed in [59j . 
particularly as a microscopic calculation [3S] having exactly the symmetry we have 
invoked reproduces the experimental data for YMn2 05 quite well. 
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